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1. Introduction. 

In recent years, the active boundary- layer control by wall-heating has 
received a great deal of attention due to its potential applications to drag 
reduction in subsonic flow over an airfoil. To support an experimental inves- 
tigation undertaken at Langley, a research program on the theoretical study 
of active flow control by surface heating/cooling has been initiated by the 
principal investigator. 

The problem of boundary-layer control has been investigated by various 
researchers for many years, (see, e.g., Chapter 14, [1]). Recently Liepmann, 
Brown and Nosenchuck [2,3] introducted the wall-heating as an active con- 
trol agent to suppress the Tollmien-Schlichting waves in the boundary-layer 
in a water tunnel. In their experiments, the boundary-layer distrubances are 
also excited by surface heating. They demostrated convincingly that, for a 
simple pluse-type of disturbances, the periodic surface-heating with properly 
adjusted phase and amplitude can effectively reduce the level of such distur- 
bances. This new technique is very promising as a tool in the boundary-layer 
control. Later on Maestrello [4] has shown that, in a wind-tunnel, localized 
surface heating can be used to trigger instantaneous transition or to reduce 
the level of disturbances by coupling it to a feedback acoustic device. To pro- 
vide a theoretical justification, Maestrello and Ting[5] performed an analysis 
of active control by surface heating. 

So far the laboratory experiments have been conducted under ideal condi- 
tions, where, for instance, the flow distrubance is known to be simple periodic. 
The control problem becomes relatively easy. To realize such active control in 
a real flow situation with unpredictable disturbances, one needs an automatic 
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control device. Therefore a rational formulation of such a problem should be 
based on the optimal control theory [6], which has been the guiding principle 
in the course of our investigation. In contrast with other analytic studies 
[7,8], we seem to be the first ones who introduced the optimal control con- 
cept into the problem of active flow control by surface heating. A complete 
formulation of this problem would involve a coupled system of Navier-Stokes 
and energy equations subject to a boundary control by surface-heating el- 
ements. However, even without the control, the solution to such problem 
is already formidable. Clearly one has to start with a simplified model and 
then approaches the original problem step by step in improving the physical 
approximation . 

Since the theory of optimal flow control is relatively new, the primary 
goal of this research project has been to explore the potential application 
of the optimal control theory, which is well developed for finite-dimensional 
systems, to a fluid dynamical system. Such system, as an example of the 
distributed-parameter system, is infinite-dimensional and, therefore, much 
more difficult to deal with. In this report, we shall present a general theory 
of the boundary Layer Control by surface heating, which is given in Section 2. 
In the next Section, we will describe some analytical results for a simplified 
model, i.e., the optimal control of temperature fluctuations in a shear flow. 
The result may provide a clue to the effectiveness of the active feedback 
control of a boundary layer flow by wall heating. In a practical situation, 
the feedback control may not be feasible from the instrumentational point of 
view. In this case the vibrational control introduced in systems science can 
provide a useful alternative. In Section 4 we briefly explain this principle 
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and applies it to the control an unstable wavepacket in a parallel shear flow. 
Application of such novel control technique to more complex fluid mechanical 
systems will be explored further. 

2. General Theory of Boundary Layer Control by Heating 

Consider the viscous flow over a semi-infinite plate lying on y = 0, x > 0, 
in the x — y plane. Let the upstream U^x) be perturbed by a disturbance 
U(x,y,t ) so that the upstream velocity components are 


u 


— Uo o + eU 0 (x,y,t),v = 0 


where t is a small parameter measuring the magintude of the disturbance. 
Suppose u 0 ,u 0 ,p 0 and 0 o are the Blasius velocity components, the pressure 
and the temperature, respecitvely. In the boundary-layer, we set u — u Q + 
eu,v = v 0 + ev,p = po + ep and 0 = 0 o + t0. Then, by a linear stability 
analysis of the non-steady thermal boundary-layer equation, we get 
, , du du du du 0 du 0 , „ * , , 

(1) It + “°& + + + ~Sy v + i 


( 2 ) 


,86 86 86 , , d 2 6 ^. 8u 0 .. 8u . 

C (~dt +u °^ - k 7^ + 2 ^ a„ 


’8x 


8y 8y 2 


8y dy 


( 3 ) 


du 8v_ 
dx^dy ’ 


with the boundary conditions 


( 4 ) 


(u,u) = { 


(fi 9)i 1 = 0 

(0,0), 7/ = 0, a: > 0 

(U 0 , 0), x < 0 or y = 00 
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(5) 


9 = 9 0 , for t - 0, and at x = 0 or oo for t > 0, 


39 

k 0 -r- = —q(x,t) at y = 0 

oy 

In (1) the kinematic viscosiy y. is assumed to depend on the temperature and 
the source term h is given by 

, x ,/ \ 1 dp 1 1 dUo i TT dUo , TJ dUoo x 

( 6 ) h(x,yJ) = --g i =--(-g r + V + 

as a result of the boundary-layer approximation (see Chapter. 15, [1]). In 
(2), c is the specific heat and the thermal diffusivity k is assumed to be 
constant. The heat flux q in (5) is the control. 

As a first approximation, we follow a reasoning by Liepmann et al [2,3]. 
Since 


(7) 




which is proportional to the heat flux across the wall and — 0 away from 
the wall, thus we have j- y {y |j) = +M0, where is an effective normal 

velocity given by v e = — (|^) = q near y = 0, and v e = 0, away from y — 0. 
Introduce the differential operator 


(8) 


,d d d 
Au = { di + Uo te +V °te~ fi ‘ )U ' 


where A = By the above physical argument, the system (1) - (5) 

can be replaced by 

du 0 


(9) 


Au + 0v = h , with /? = 


dy' 


(10) 


3n 3v 

— — f- — = 0 for x > 0, y > 0 and i > 0, 
ox oy 
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K«) = { 


{f,g ) at i = 0, 

(0,u e ) at ?/ = 0, a; > 0 
(£/ 0 , 0) for x < 0 or y = oo 


where the energy equation (2) may be neglected. 

The above systems form the state equations for the control problem where 
the effective velocity v e = q will be used as a boundary control. 

Let D = {x > 0, 7/ > 0} be the region above the plate; T 0 and I\ be 
two disjoint parts of the wall T = {x > 0,y = 0}. For a feedback control, 
we measure the shear-stress distribution t over Tj to yield the observation 
equation 

. du 

(11) T(x,t) = pn— on lx. 

For the optimality criterion, we introduce an objective functional or perfor- 
mance index J(v e ). For example, one may choose 

(12) J(i> e )= [ T [ T 2 (x,t)p(x,t)clxdt + [ T [ v e {x,t)-{(x,t)dxdt : 

Jo jTi Jro 

where the second integral represents the control cost function, while p and 7 
are appropriate weight functions. 

Let Q be the admissible set of controls v e of functions of t and x , for 
0 < t < T and 1 in To (the porous region). Then the optimal control 
problem can be stated as follows: 

Find an optimal surface- heating rule v*(r) from the admissible class Q 
such that v* minimizes the observed wall-shear and the control cost, i.e. 

(13) J«) = min J(v e ). 

v * m Q 
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This can be formulated as a problem in variational calculus. If such a solution 
q* can be constructed, it will yield a feedback control law 

(14) v e = <(t). 

This control law will regulate the injection velocity distribution automatically 
based on the wall-shear input. 

To be specific, let us consider two special types of problems: 

1). Feedback Control of Normal - Mode Instabilities. 

For a parallel shear flow subject only to an initial perturbation, we set 
v 0 = 0 and h = 0 in Equation (8). By introducing a stream function ip, the 
system (9) - (10) yields 

(15) \j t + Ml = "o (!/)£<? + R~' A V. 

p(x,y,0) = <p 0 (x,y), 

r 

|y=0= Ve( x > 0) ajj'r’ lj/=0 = 0, 

( 16 ) j 

= 0 as y — ► oo, 

where R is the Reynolds number. 

Given the observed wall-shear r, the feedback control law is assumed to 
be of the form: 

(17) v e (x,t) = J g(x - t)r({,t)dt = p J g(x - 

where g is an optimal transfer function as 3 ), et to be determined. For a normal 
mode analysis, let 

(18) <p(y, t, k) = J v?(x, y, t)e' kx dx = 8(y, k)e ,wt 
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Then the system (15) - (16) reduces to an eigenvalue problem for the Orr- 
Sommerfeld equation with a modified boundary condition: 


(19) LO = XM0,y > 0 

ikO = ggO" and O' = 0 at y = 0, 
0 = O' = 0 at y = oo. 


Here we have put 

L = -(DW)\D = £, 

M = (U - c)(D 2 - k 2 ) - U\ U = U 0 /c] c = w/k, 

A = ( — iw), 

and 

g(k) = J g(x)e ,kx dx. 

Since the real part Re {A} of the eigenvalue A yields the growth rate (if 
positive) for the normal mode with wave number k. Therefore we choose the 
performance index J(v e ) = Re (A(u)}. The optimal control problem is to 
determine the transfer function g (or g) so as to minimize the growth rate 
J(v e ) = J(g). The optimal transfer functions g* will give the optimal control 
law via Eq. (17). 

2). Optimal Control of Externally Excited Instabilities. 

In constrast with the previous problem, the external perturbation h in Eq. 
(8) is non- zero. By method of superposition, one assumes that the normal 
mode (18) and 

(20) h(y,t,k) = J h(x,y,t)e ,kx dx = h(y,k)e~ twt . 
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In view of (20), a Fourier transform of the system (15)-(16) yields a nonho- 
mogeneous Orr-Sommerfeld equation with a modified boundary condition: 

(21) (. L-\M)0 = h{y,k),y>0 , 
ik9 = ggO" and O' = 0 at y = 0, 

0 = O' — 0 at y = oo, 

which is a nonhomogeneous boundary- value problem. For a feedback control 
of the form (17), we choose the mean kinetic energy of the pertubation as 
the performance index: 

J{v e ) = Jirn ~ J y, t) + v 2 y (x, y, t)}dt dx dy, 

which in terms of 9, reduces to 

/•OO _ 

( 22 ) j(v e ) = I { I e Xv, k ) I 2 + fc2 1 %»*) I ) d y ■ 

J 0 

Since J(v e ) = J(g), the control problem at hand is to determine the transfer 
function g* which minimizes the mean kinetic energy J . This problem may 
be solved by the method of Green’s function. Let G(y,ri,k ) be the Green s 
function for the boundary- value problem (21) so that 

fOO 

(23) 0(y, k) = G{y,r},k)h(ri,k)dri. 

Note that, through a boundary condition in (21), G depends smoothly on 
g (or g). A substitution of (23) into (22) shows that J{g ) is a functional 
of g. Therefore the optimal g* (or g*) can be determined by the variational 
equation 

(24) SJ(g) = 0. 
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Thereby the optimal control law can be determined. For this problem the 
major task is to carry out the variational analysis. Then numerical compu- 
tation may be done to show the effectiveness of the optimal control and its 
dependence on various physical parameters. 

3. Solution of a Simplified Control Problem 

As mentioned before, due to the complexity of the surface heating control 
problem, an analytical solution is unattainable. Even the numerical solution 
is difficult. In this section we shall present a simplified model problem which 


can be analyzed rather completely. 

Recall that the analysis of heated boundary layers is based on a coupled 
system of the momentum equation (1) and the energy equation (2). One 
notes that the coupling between these equation is due to the dependence of 
the viscosity on the temperature. It seems safe to say that, in the absence 
of mechanical disturbance, the flow instability can be achieved by reducing 
the temperature fluctuation, which is the only source of external excitations. 
Furthermore we regard the dynamical heat production term in the RHS of 
Eq. (2) as an external perturbation denoted by, say, cf(x,y,t ), and assume 
that the flow is parallel to the plate with «o = U(y) and vq = 0. Then we 
are led to considering the following optimal heat regulating problem: 




d 2 d 2 \ 

w + di *) 9 + f ( x,y ' t ^ 


(26) k 0 -(x,o,i) = -q(x,t), 

dy 

for t > 0, — oo < x < oo, and y > 0, where v = k/c is the thermal diffusivity 
and / is an equivalent source of thermal disturbance. Of particular physical 
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interest is the case where the perturbation is persistent in time and localized 
in space. Thus the perturbing source function / is assumed to satisfy the 
mean-square integrability condition 

1 ft fOO fOO 

(27) lim - / / / \f(x,y,t)\ 2 dt dx dy < °o. 

T —>oo ± J o J-oo JO 

As a consequence, the heat flux q , the active control, is expected to have a 
similar property 

1 fT foo 

(28) lim - / / \q{x,t)\ 2 dt dx < oo. 

T— ► OO i Jo J — oo 

As an optimality criterion, the objective function J given below will be min- 
imized, 

1 fT foo „ 

(29) J{q) = lim — f / {|0(a;,(M)| + N\q(x,t)\ dt dx 

T—*oo 1 JQ J-oo 

where the first term in the integrand yields the mean-square temperature 
fluctuation along the wall, while the second one is a measure of the mean- 
square control cost, where iV > 0 is a cost parameter. For N 0, the 
minimization of J gives the optimal control q which is the best possible for 
the system to realize. On the other hand, for > 0, the optimal solution 
q will be to minimize the wall temperature fluctuation without excessive use 
of the control action. The choice of the number N depends subjectively on 
the relative importance assigned to two competing factors in the optimality 
criterion. Now the mathematical problem for the control of thermal distur- 
bances can be formulated as follows: For a given heat source perturbation 
/(x,j/,i) satisfying the condition (27), find the optimal control with 

the property (28), which minimizes the objective function (29) . Here the 
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“state” 0 is the solution of the system (25) and (26) subject to the zero initial 
condition. Note that, since the transient part of the solution to the system 
will be wiped out by the time average in (29), without loss of generality the 
initial condition may be taken to be zero. 

Since the governing equation is linear and the functional J is quadratic, 
the problem can be solved by the principle of superposition. Therefore it is 
sufficient to consider a time-harmonic perturbation 


(30) 


= g(x,y)e ,wt . 


Then it is possible to seek a time-harmonic solution 
(31) 9(x,y,t) = 4>(x,y)e~' wt , 


and 


(32) 


q(x,t) = r(x)e 


—i uft 


A substitution of (30)-(32) into (25), (26), and (29) yields 


(33) 


L<j> = 1/V 2 0 - U(y)^ + iu<t> = - g(x,y ) 
ax 


(34) 


k 0 ^(x,o) = ~r(x), 

oy 


/ oo _ 

{|^,o)| 2 + iV|r(a:)| }dx, 

-OO 

where <j> vanishes as |x| — > oo or as y — * oo. It is well known that the 
necessary condition for J being minimal is 


(36) 


6J{r) = 0, 
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where 8 J means the variation of J with respect to r. Note that the quantities 
(j),r are complex. Denote their complex conjugates by <f> * and r* respectively. 
By introducing an adjoint state ?/> to <j>, one can derive the following optimal- 
ity system [6] 

(37) L(j> = uV 2 (f> — U(y) ^ + iu<f> == -g{x,y), 

(38) L*j> = + U(y) ^ - iwtf = 0, 


(39) 


, d<f> l . . 

k °~ dy ^ x '° 2 ~ = °» 


(40) <f>(x,o) + k Q ^(x,o) = 0. 

It is seen that, in order to find the optimal control r, one must solve an 
extended system of coupled equations (37)-(40). Then the optimal solution 
is given by 


(41) 


r - r(x) = ~—ip(x,o) 


and 


(42) <j(x,f) = ?(x) e xwt 

gives the optimal control for the time-harmonic disturbance (30). In the 
multiple frequency case, the optimal control q is a sum of the single-frequency 
solutions by superposition. 
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Clearly, by a Fourier transform in x, the optimality system (37)-(40) is 
reducible to a one-dimensional problem: 

(43) v<j> n — [i^A 2 + i{XU T = J/)j 

(44) vi>" - [ 1 /A 2 - i(XU + u))]i> = 0, 

(45) kj(\,o)-j r i>(\,o) = 0, 

(46) k 0 tp'(\, o) + <£(A, o) = 0, 

where the Fourier transform /(A) of a function /(x) is defined as 

/(A) = / f(x)e lXx dx. 

J — CO 

Let 4> p be a particular solution and two linearly independent comple- 
mentary solutions of the equation (43) such that \<j> 2 \ * 00 as V * Tb en 

the complex conjugates <j>\ and are linearly independent solutions of (44). 
Therefore the bounded solutions of (43) and (44) are given by 

(47) #(A,y) = ^ P (A,y)-a^x(A,y), 

(48) $(A,y) = ,y), 

where, in view of the boundary conditions (45) and (46), the constants a and 
b are found to be 

a(X) = + jji*,i p }(X,o)/d(\), 

6(A) = k a (M' p -^i r )(\,o)/d(X), 
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where 


4A) = (^;| I + ^|*|)(A,o). 

In view of the equations (43) and (48), by an inverse Fourier transform, the 
optimal control law is obtained, 

(49) ^( X ) = _ 27t nI K^d^°) e ~ xXx ^- 

Under the optimal control, the state of thermal fluctuation is given by the 
inverse transform of (47): 

1 f°° 

(50) <j>(x,y) = — e~ lXc [<j>p{\,rj) - a(\)4>i(\,y)]d\. 

ZlT J — od 

On the other hand, when there is no control (r = 0), the temperature fluc- 
tuation <j>o{x,y ) is given by 

1 rOO ^ a. 

(51) <f> 0 (x,y) = —[ e~ ,Xx {<t>p{\, y) - a 0 (A)^i(A, y)]d\, 

ZlT J-OO 


where 

a oW = 

To assess the effectiveness of the control, we may introduce either the uniform 
reduction ratio 


(52) 


e 0 = maxe(x) 


max 


4>{x, o) 

4>o(x,o) 


or the (root) mean reduction ratio 

( roo roo ^ 1/2 

(53) £ m = [J J<f>{x,o)\ 2 dx I j_JMx,0)\ d*j 

( roo a r°° a r, 1 */ 2 

= [J J^x,o)\ 2 d \/ , 
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by Parseval’s equality in Fourier transform. To give a qualitative physical 
interpretation of our results, it is instructive to go through some examples in 
detail. 

Two examples corresponding to a time-harmonic perturbation f(x,y,t) 
of the form (30) will be considered. In particular the spatial distribution 
function g is assumed to be 

(54) g(x , y) = A exp{i/ca: - ^ai 2 - /3y} 

which represents a decaying surface-wave disturbance with amplitude A , the 
wave number k, and the exponential decay parameters a and j3. Then the 
Fourier transform of g in x is given by 


(55) 


<K A >2/) = <7o(A)e Py , 


where 

(56) 


9o 


(A ) = A\K-^ x+ ^. 
a 


As the first example, let the velocity U = U 0 be a constant. Then the 
transformed optimality system (43)-(46) with g given by (54) can be solved 
analytically. The transformed optimal solutions (47) and (48) become 


(57) 

(58) 
where 

(59) 




e~ 0y - 


T + Nkl/3? ' 


0 -t*y 


#(A ,y) = h(\) 


Nk 0 {{ - 13) 


Li + ^lkl 2 ] 


„-(*v 


1 (\\- 

CW-P 2 ' 
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and £*(A) is the complex conjugate of £(A) defined by 
(60) £(A) = {A 2 + i(\U a + u/v)Y'\ Re{({ A)} > 0. 


By the inverse Fourier transform (49), the optimal control law is found to be 

, K r gJX )e-“‘ 

(t.1) r x — 2*y-»(l + JVi*|«A)| 1 )K(A) + /Jl’ 

Along the wall, from (50) and (51), the controlled and the uncontrolled tem- 


perature fluctuations are given by 

Nkl r°° £*0)&>( A) 


r< 

(62) <f>(x,o) = ^J_ 


2 tt j-ooitw+m+Nwam 


l~ iXx dx , 


(63) 


Mx ' o) = bL 


9o{ A) 


— t A x 


d\. 


«A)K(A) + jS] 

One observes that, if N is large, then ;'(x) ~ o and <^(x,o) ~ <f> 0 (x,o). That 
is, when the control is costly, no control action is taken so that the state is 
unchanged. While, if the control cost parameter N is small, we get 

9o( A)e" ,Ar 


(64) 


'(*) ~ J_ 


CO £(A) + /? 


■dA, 


and <f>(x,o) ~ 0. Thus the unrestricted thermal control almost completely 
eliminates the temperature fluctuation along the wall. In fact, in this case, 
the control action is to cancel out the heat flux across the wall at any cost. 

To simplify the results, let the decay factor a be small. In view of (56), 
as a 0,j approaches a <5-function and the otpimal control (55) reduces to 


(65) 


f(x) ~ — 


Ak 0 e ,KX 


({(-*) + 0(i + ' V *2I«-*)I 2 ) 


16 


The corresponding residual temperature fluctuation (51) becomes 


( 66 ) 


ANk 

* {x,o) ~ [«-*) + m + /w-» 2 I{(-k)I 2 )' 


Recall, by noting (54), that 


(67) i{k) = {k 2 + i(kU 0 + u/v} l/2 , Re{() > 0, 


which gives 


|i(fc )| 2 = {A-" + (i-t/ a + us) 1 /^}*/ 2 . 


As it turns out, in this case, the uniform reduction ratio and the mean 
reduction ratio, defined by (52) and (53) respectively, coincide. In fact we 
have 


( 68 ) 


Nki\g-K )\ 2 

i + Nk 2 \a -*)\ 2 


< i, 


where £ = £ 0 = £ m . 

From (54) one sees that, as a — + 0, the excited disturbance is a surface 
plane wave. The otpimal response t in (65) is also of the same wave form with 
a phase-shift arg.{£(-K) + /?}. Under the control, the surface temperature 
fluctuation is reduced to (66), which has the reduction factor e given by (68). 
It shows clearly that £ increases from 0 to 1 as tj = Nkl\£(— k)| 2 increases 
from 0 to oo. For a fixed control cost parameter N, rj is an increasing function 
of k and u;, but a decreasing function of v. Therefore the reduction ratio £ 
and, hence, the effectiveness of the thermal control decreases as the wave 
number k or the frequency u of the distrubance increase, while they increase 
as the thermal diffusivity v of the fluid increases. These results are physically 
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plausible and may provide some insights into a certain flow control problem 
by surface heating or cooling. 

As the second example, assume the shear velocity U varying slowly from 
zero at the wall to the free stream velocity U 0 at the height y = |, where 8 is 
small. Then, by applying the multiple-scale or the ^VKB Method [9], one can 
construct asymptotic solutions to the optimality system (43)-(46), by using 
8 as the small parameter. Without giving the derivation, it can be shown 
that the first-order asymptotic approximation yields the complementary and 
particular solutions to equation (43) as follows: 

(69) ha(Ky) ~a{\,y)e* Q{X ' y) 
and 

(70) 

where 

(71) 

(72) 
and 

ft(A,y)= / ri{\,s)ds. 

Jo 

For a small decay factor, a ~ 0, by some tedious by straight-forward al- 
gebraic manipulations, one can obtain the otpimal control f and the residual 


<^ P (A, y) ~ [2t/(A,o)] x {^i(A ,y) J q h{^,s)g{\,s)d& 
+ <fo(A, y) J <^i(A, s)<7(A ,s)ds} 


ri(\,y) = {A 2 + i v~ l [\U{y) + cj ]} 1/r2 , Rey > 0, 


1/2 


a(x ’ y) = teJ ’ 
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wall temperature fluctuation <f>(x,o ) in a simple form: 


(73) 


r(x) ~ —A 


(i + NkSM-K.o)! 2 )' 


(74) 

where 


<f>{x,o) ~ 7 (-k) 


Ic 2 0 Nti*(-k,o) 


1 + k*N\ti(-K,o)\* 


e ,KX , 


(75) 


7(A) = f°° a(\,T)e- Q{X ' T) - 0T dr 
J 0 


From (66) it follows that the reduction ratio 

™ . _ Nkl\ri{-K,o) j 2 

1 j l + iV^|, 7 (-/c )0 )|2- 

In contrast with the results (65)-(6S) for a uniform shear profile, the 
above results (72)-(75) show that, for a slowly varying shear flow, the velocity 
variation introduces a shape factor 7 as defined by (69) which modifies the 
control and the corresponding temperature fluctuation. In fact it is easy to 
verify that if it = U 0 , the results (72)-(75) reduce to the previous ones (65)- 
(68) with a uniform shear velocity, as they should. Therefore the previous 
physical interpretation of those results is still valid, at least qualitatively. 

To issustrate the results graphically, some numerical calculations have 
been performed when the shear velocity profile is uniform. The numerical 
results are displayed in Figs. 2-8. Figure 2 shows the periodic variation 
of the thermal control input along the wall as the nondimensionalized wave 
number k varies from 0 to 2, where 77 = i?e{r}, and A = k 0 = f3 = N = \. 
It is interesting to note that, in response to the periodic distrubance, the 
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controlling heat input rises steeply, as k increases, and then falls off sharply. 
This trend is clearly seen in Fig. 3. In Fig. 4, two sets of curves correspond to 
the controlled wall-temperature distribution (Ci = Re() and the uncontrolled 
one (£1 = Re(^ 0 ), for different values of k. The associated optimal control is 
given in Fig. 2. The dependence of the amplitudes of the controlled and the 
uncontrolled wall-temperature fluctuations on the wave number k is shown 
in Fig. 5. The above two figures show clearly the effect of the thermal 
control on the wall temperature fluctuation at a moderate control cost factor 
a = Nk 2 0 = 1. With k = 1 and the rest of the parameters fixed as before, the 
set of curves in Fig. 6 displays the residual wall-temperature distributions 
under the optimal control for several values of a, while the corresponding 
amplitudes are plotted in Fig. 7 against the effective control cost parameter 
a for several values of k. Finally in Fig. 8, the reduction ratio £, which 
measures the effectiveness of the control, is sketched as function of the wave 
number «, for several values of a. Obviously the optimal control is most 
effective in reducing the thermal disturbances with small wave numbers and 
a low control cost parameter. The control loses its effectiveness for short-wave 
disturbances or when the control cost parameter gets too high. These findings 
seem to be consistent with one’s physical intuition. Also it is interesting to 
note that, for small wave numbers, the reduction ratio dips into a minimum 
value before it takes off and increases monotonically to one. 

4. Vibrational Control of Unstable Wavepackets 

The principle of vibrational control of dynamical systems was proposed by 
Meerkov [10] for control problems where feedback and feedforward principles 
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vibrational control ^p(wt ) with a high frequency u? >> 1. To this end, let us 
scale the time by setting s = wt and e = 1/cj. Then Eq. (77) can be written 
as 

(79) ^ = eF{x,t,<p{s)), 

where t = es becomes a slow variable. Note that Eq. (79) is in the canonical 
form, in the language of nonlinear oscillation theory. Thus, by the averaging 
principle of Bogoliubov and Mitropolski [13], the average state y{t) satisfies 
the equation: 

(80) = F(y,t), 
where 

F{y,t)= [ F(y,t,<p(s))ds 
Jo 

is the average of F is s over a period with y,t held fixed. By assumption, 
y = x o is also an equilibrium point of Eq. (79). If it is possible to find a mean- 
zero periodic excitation at a high frequency u in Eq. (78) so that the 

average equation (79) is stable at xq in the mean, then the system (77) with 
y, = 0 is said to be vibrationally controllable near x = zo- There are two ways 
to excite the system: additive excitation and multiplicative (or parametric) 
excitation. In the following application, we shall be concerned with the case 
of additive vibrational control only. In this case F(x,‘,y) = Fo(x,t) + y so 
that Eq. (77) reads 

(81) ^ = F 0 (x, t) + y. 
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For vibrational control, we set f.i = in the above equation 

to get the controlled system: 

dx 

(82) — = F 0 {x,t) + <p{ut). 

Now let 

x(t) = (p{ut) + C(0» 

which is then substituted into Eq. (82) to get 

(83) = Fo(<p(ut) + C,t). 

If we set F 0 {<p(ut) + (,<) = then Eq. (83) becomes a special 

case of the system (78) and the vibrational control method described above 
in applicable to this problem. In fact, since has a zero mean, the 

controlability of Eq. (82) near x 0 implies the same for Eq. (81). Now we are 
going to apply this technique to control the instability of wavepackets in a 
shear flow as announced. 

Consider a parallel shear flow in the x-direction. For a two-dimensional 
problem, let u(x,t) be the velocity perturbation about the mean flow U at 
the Reynolds number R near the critical value R c . Let us write 

R = R c ± 6 2 , 8 = 1 R — R c |^ 2 

and 

u(x,t) = 2Re{A((,T)e i ( a '*-*‘ t) } 

where £ = 8x and r = 8 2 t , that is, we are seeking a solution in the form of 
slowly varying wavepacket. It was shown that the variation of the amplitude 
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A of the most unstable mode satisfies the equation [12]. 


x dA dA . . 1 c2/>\ a a a 

( 84 ) Ji+ C ’d( 2 ~ i SHW A ' 

where c g is the group velocity; i is the Landau constant, and the complex 
constants c and a appear in the following expansion: 

A = —iu = iu> c — ic g (a — a c ) ± kS 1 ^ 2 — a(a — a c ) 2 + 

To simplify Eq. (84), we define 

t 7 = (f — c 3 t ) and B = A/ 6, 


so that it yields: 
(85) 


a i- a w = kB ~>' 1B - 


It is known that a and k have a positive real part and that, if R > R c and 
l r = Re {£} < 0, the amplitude may tend to infinity in finite time. Such a 
burst is reminiscent of sudden transition to turbulence, (see p 452, [12]). In 
order to control the instability, we introduce an additive vibriational control 
1/7(77, wt ) and consider the controlled equation: 


(86) 


^ - a = kB - \m 2 B + VO lM)- 


dr drjt 

Let us consider a plane wave disturbance of the form 


B{ 77, r) = b(T)e ,/3r> , Im /? = 0, 


and the corresponding control 


1/7(77, ut ) = . 
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Then Eq. (86) yields 

(87) ^ = {k-a/3 2 )b-^£\b\ 2 b + ^>(ujT), 

which is of the form (82). If there is no control (y> = 0), we can derive from 

(87) the following equation for R = |6| 2 : 

(88) ^ = 2<j R- l r R\ 

where a = Re {& — a/3 2 ) and l r — Re {^}. The logistic equation (88) has the 
solution 

(89) R{t) = cR x e 2aT j{l +ce 2aT ), 

where c = , R 0 = R{ 0) and R x = From (89) it is clear that, if 

Re a > 0 and R\ > Rq, the equilibrium R = 0 is unstable. On the other 
hand, if Re a < 0 and R x > R 0 , R = 0 is asymptotically stable. However, 
for any cr, a burst may occur if R\ < Ro . Now returning to the controlled 
equation (87), we let 

(90) 6(r) = <p(ut) + b £ (s), 
where s = t / e, e = - . Then b c satisfies 

ft CJ 

% = + <■')}. 
which yields the following average equation: 

(91) 

where (p 2 = / 0 1 ^ 2 (^)c ?5 and (p 3 = 0. In conti'ast with Eq. (87) with p = 0, 
one cannot solve for R = \b\ 2 in a closed form. To see the effect of the 
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vibrational control, let us consider the linear stability near 6 — 0. Writing 
6 r = Re 6, 6, = Im 6, etc., the linearized equation of (91) can be written as: 



where 

7l = (°V - ^r^ 2 ),7l(^r - 

and 

72 = (<7. - |Av , 2 ) > 72 = (o-i - h.<£ 2 )- 

Let us substitute for (6 r ,6^) by (a r ,a,)e AT in Eq. (92) to yield an eigenvalue 
problem, for which the characteristic equation reads 

(A — 71) 72 
-72 (A — Tl) 

The roots of the equation are 

A, 2 = (a, - l r tf) ± - (ffi - <,V 2 ) 2 } 1/2 . 

By the method of finding a maximum in calculus, it is easy to show that 

Re A < {(J r - i T p 2 ) + {\\ItP? + ^, 2 } 1/2 - 

Thus, if c i — 0, the equilibrium 6 = 0 is stable provided that a T < \i T (p 2 ■ 
In contrast with the stability condition a r < 0 for the uncontrolled case, the 
vibrational conti’ol can stabilize the system if we make <p 2 large enough. For 
any <t,(^ 0), this is also true. Of course we pay a price for stability by allowing 
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* 


a small zero-mean oscillation about the equilibrium. Further applications of 
vibrational control to fluid mechanical systems will be discussed in our future 
work. 
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